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We propose a complete quantum computation process on a topologically protected Josephson 
junction array system, originally proposed by Ioffe and Feigel'man [Phys. Rev. B 66, 224503 
(2002)]. Logical qubits for computation are encoded in the punctured array. The number of qubits 
is determined by the number of holes. The topological degeneracy is lightly shifted by tuning the flux 
along specific paths. We show how to perform both single-qubit and basic quantum-gate operations 
in this system, especially the controlled-NOT (CNOT) gate. 
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I. INTRODUCTION 

Quantum algorithm shows us a tempting future to con- 
quer some problems that are hard to solve using classi- 
cal computers [lj, such as, factorizing large numbers [2j, 
searching in a large database[3], and so on. However, the 
physical realization of a quantum computer turns out to 
be a really difficult work. The greatest obstacle is the 
decoherence of system induced by unavoidable coupling 
with the surrounding environment Since quantum de- 
vices work in a microscopic scale, their interactions with 
the environment become fatal to the quantum coherence 
of the system. Someting has to be done to increase the 
decoherence time so as to complete enough operations 
within so short a peroid of time. 

Therefore, error correcting code (ECC) is neededp]]. 
With the help of ECC a small error rate can be tolerated. 
Howerver, since the necessary error-rate threshold still 
lies far beyond the present available technology, other 
methods to fight against error should be considered. One 
promising method is topological quantum computation 
(TQC) 110. 

Topological order[8] [9] [10] [1 1] is a new kind of quantum 
order. Although the mechanism of topological order is 
not quite clear, we can find the following features in a 
system in topological order: 

1 . The degeneracy of the ground state closely relates 
to the topology of the physical system (usually 2- 
dimensional) . 

2. The system is immune to local perturbation and 
this property can be used to protect qubits. 

3. Excited quasiparticles may obey anyon statistics 
(Abelian or non- Abelian) Q . 

With its great advantages, non- Abelian TQC[7] attracts 
us the most. In a non- Abelian topological phase, not 
only can the encoded qubits be protected, but also the 
operation process can be protected intrinsically through 



the braiding of non- Abelian anyons. Much has been dis- 
cussed about the Fibonacci anyons that can be used for 
TQC directly in Quantum Hall Effect. However it is hard 
to control anyons in real quantum Hall systems. Be- 
sides, Kitaev proposed a spin-1/2 honeycomb model[T^|. 
which seems achievable, followed by some different pos- 
sible physical implementations[I3] [14] [1 5| [l6l| [l7j). The 
work in Ref.flT} proposed a more sophisticated method 
for the anyon detection process. Unforturnately, the non- 
Abelian anyons that emerged in this model do not act 
well enough to complete the full braiding operation of 
TQC. Attractively, B. Doucot, et al. gave a complete de- 
scription of the rhombus Josephson junction array (JJA) 
with Abelian and non- Abelian gauge theories [la ]. How- 
ever, this is too complex for implementation. 

Because of all these difficulties we are facing, we will 
turn to the Abelian topological phase in which qubits 
can be protected as well. Recently, Kou showed how to 
do computations in the toric code model by using the 
tunneling effect theoretically [l9j [20j |2l| . However, the 
original toric code model calls for four-body interaction, 
which is hard to implement in real systems. 

The Josephson junction seems to one of the candi- 
dates to implement a realistic quantum computation 
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24] [25]. Ioffe and Feigel'man proposed a hexago- 
nal Josephson junction lattice that showed some topolog- 
ical protecting properties which can be used as quantum 
memory [2oT |. 

Based on their work, we propose a complete TQC 
scheme in this article. We focus on encoding logical 
qubits in the array and the construction of quantum 
gates. Similar to the method of Kou, we also show how 
to implement fundamental quantum gates, especially the 
CNOT gate. 

The article is organized as follows. In Sec. II, we intro- 
duce Ioffe and Feigel'man's model in detail. In Sec. Ill, 
we study the ground-state behavior. Topological degen- 
eracy is used to encode qubits in the array. In Sec. IV 
we show how to manipulate the qubits to complete basic 
quantum gate operation. Finally, we draw summary in 
Sec. V. 
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Figure 1: (Color online) Construction of the basic element in 
the array. Four identical Josephson junctions form a closed 
loop, embedding flux $o/2. The classical potential is shown 
on the right. 



II. HEXAGONAL JOSEPHSON LATTICE 

To describe the computation process we first introduce 
the model proposed by Ioffe and Feigel'man [25| in this 
section. In this work, They showed that the array system 
exhibits some topological protecting properties and thus 
can be used as quantum memory. We discuss how to do 
computations based on their work in the next section. 
Now we recount some basic properties of their model in 
detail first. These properties are essential for the con- 
struction of computation gates. For the convenience of 
the computation we talk about, we make a small modifi- 
cation with the boundary and discuss further the quasi- 
particle excitation, especially the vortex, which agrees, 
in essence, with the original model. 

First, we give the basic building element and explain 
how it works. Then we show how the array is arranged 
and give the dynamics of the system. Finally, we talk 
about the topological degeneracy of the ground space and 
the low-lying excitations. 



A. Basic rhombus 

As shown in FigJTJ the basic building block is a rhom- 
bus made of a four-junction superconducting loop placed 
in a magnetic field. The flux through the rhombus is 
$o/2, where $0 = h/2e is the flux quanta. We re- 
quire that all the Josephson junctions are the same and 
work in the "phase regime", that is, Ej Ec, where 
Ej = (h/2e)I c and E c = e 2 /2C. 

This system, as is much discussed in the literature, 
forms a simple two-level "flux qubit'' (or persistent- 
current qubit) [23[ [2(| . In the following, we briefly explain 
how the qubit system emerges. 

We denote the phase of the superconductor islands 
with ipi. The Josephson system is described with the 
Lagrangian 



Here, we make an approximation that three of the junc- 
tions have the same phase difference. As all phase differ- 
ences should sum up to 27r$ ea;t /$0j we get 
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where A<fi indicates the phase difference between the 
two far ends of the rhombus. Neglecting the kinetic 
part, this energy has two equivalent classical minima at 
A<p = ±7r/2 (as shown in FigJT]), that is, the phase dif- 
ference across each junction is ±7r/4. So the two clas- 
sical ground states correspond to states in which there 
are both clockwise and counterclockwise supercurrents in 
the loop. This is the reason why we call it a persistent- 
current qubit. We denote the two states as | f) and | J,). 
Quantum fluctuation makes the electron tunnel between 
the two classical states in the wells and that splits the 
degenerate ground states. The low-energy configuration 
can be regarded as a two-level qubit system. The effec- 
tive Hamiltonian is 



H = ta x , 



(3) 



where t ~ ^/EjEce~ s ° , and Sq is the classical action 
taking the path between the two minima. The value of 
tunneling energy is estimated regarding the potential as 
quadratic. 

Let us see what happens when the condition is slightly 
changed from the ideal. Notice that the dynamics require 
that the flux through the rhombus is exactly half a flux 
quanta. If the external flux is a bit different from the 
ideal value, <&o/2 + 8&, we add a first order change to the 
potential, 
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As a result of the slight change of the flux, the two min- 
ima of the double well are no longer degenerate but sepa- 
rated by a gap AE = 2^/2ttEjS^/^o. Then the effective 
Hamiltonian is H = ta x + ^AEa z . 

This is important for operations. We can control the 
evolution of the system by adjusting the flux through the 
rhombus. Below, in the Josephson array building from 
the basic rhombi, logical qubits are encoded for computa- 
tion. Some gate operations are done just by the method 
of tuning the flux of the rhombi along specific paths. 
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(ipi - ipj) 2 + Ej cos(V>; - ipj - aij), (1) 



where ay = ^ f. A - dl. The quantity in the cosine term 
is the the gauge-invariant phase difference across each 
junction. 



B. Array 

Now we have got the bricks we need. In this part we 
show how to build the Josephson array system. 

The basic rhombi are arranged as shown in FigJ5] For 
simplicity, a rhombus is indicated by a short line in the 
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Figure 2: (Color online) Construction of the array with the 
basic rhombus. Six rhombi form a hexagon. Each hexagram 
embeds a flux of $ s = (n s + |)$o- The array is punctured 
with several holes. The number of the holes relates to the 
degeneracy of the ground space. The boundary, as shown 
in the lower figure and indicated by a (blue) solid line, is a 
whole superconductor wire connecting all the rhombi along 
the boundary. The flux through the polygon loops on the 
boundary are $1 = ni$o and $2 = (112 + § )$o- 

figure. Six rhombi form a hexagon. The array is placed 
in a uniform magnetic field. We should carefully tune 
the shape of the rhombi, so that the flux through the 
hexagram in the center is a half- integer multiple of <f>o, 
<f> s = (n s + l/2)$o- The centers of the hexagons are 
denoted with a, b, c, and so on and individual rhombi be- 
tween the corresponding two centers with (a6),(6c) .... 
The array contains K holes. Soon we will show that, in 
fact, each hole can be treated as a single logical qubit. 

In the array, the rhombi form loops, so we should check 
the gauge-invariant condition over again. To satisfy the 
constraint condition for all the closed loops in the array, 
the sum of the phase differences around each hexagon 
should equal 27r$ s /$ [i-e., (2n s + 1)%]. This is well 
consistent with the configuration that all the A<f> a b of the 
six rhombi on the edges of a hexagon take the ground- 
state value 7r/2 (or all take — tt/2). The flip of a single 
rhombus is not allowed because it changes the sum of 
the phase differences around the loop by it and violates 
the constraint. However, we can flip rhombi of an even 
number in a hexagon at the same time (e.g., two, four or 
six). Choosing the two ground states as a basis, we can 
write the constraint in the operator form, 

Pa = U<b = 1 (5) 
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The product takes the rhombi on the hexagon loops. 

Because of the existence of the constraint, the rhombi 
embedded in the array can no longer behave as in Eq.Q, 
which does not commute with the constraint. The sim- 
plest dynamics compatible with constraint contain flips 
of three adjacent rhombi, Q a bc — IJ ab IJ bc <J ca^ sucn that 
two rhombi flip at the same time in the corresponding 
three hexagons. This gives the main contribution to the 
Hamiltonian of the system 

Ho = -r^Qabc- (6) 

abc 

Here, r ~ \fEjEc exp(— 3So) is estimated to be analo- 
gous to t in the last section. The action in the exponential 
term is three times of that of a single rhombus because 
of the simultaneous flip of the three rhombi. 

C. Boundary 

Now we talk about the boundary. We should treat it 
carefully. The boundary here is arranged a bit differently 
from Ioffe's original model. Different boundary condi- 
tions of the toric code have been discussed in detail (27|. 
In the system considered here, similarly, the boundary is 
arranged more like the "rough edge". 

Here, we arrange the boundary rhombi in such a way 
that they behave like those in the bulk, as described by 
Eq.©. Each rhombus on the boundary still stays close 
to the other two. A superconductor line connects the 
unattached rhombi at the outer and inner boundary (like 
the blue line in Fig|3]). 

The geometry of the polygon loops on the boundary 
should be carefully treated. A boundary loop has m 
rhombi, usually three or four. We want to maintain the 
dynamics and constraint of the form as before, that is, 
two rhombi must flip at one time. So the flux embedded 
should be consistent with the sum of the phase differ- 
ences of the m rhombi, i.e., (rib + t)^0j where rij is an 
integer. Then the Hamiltonian Eq.© is still kept for the 
rhombi on the boundary. 

We complete the construction of the Josephson array 
and give the low-energy dynamics. In the following sec- 
tion, we will discuss the properties of the ground-state 
space and excitations in detail. 

D. Topological degeneracy, charge and vortex 
excitation 

In this section, we discuss the properties of the ground 
space and excitations. The topological degeneracy of the 
ground space can be used to encode logical qubits for 
computation. We show two kinds of quasiparticles in the 
system, charge and vortex. We can operate qubits by 
controlling their motion in the array. 

First, let us find the ground state of the system. Ignor- 
ing the constraint, we can easily get the ground state of 
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Hamiltonian @ as \0) =\[\+) = t) + I l»/V2. Of 
course, it is illegal and incompatible with the constraint. 
We can project it into the physical space with the help 
of an operator 



1 + Pg 

V2 



(7) 



\g) = V\0) is a ground state of the system, but not the 
only one. 

To find all the ground states, we can turn to the help 
of a path operator 



(7*) 



ib '■ 



(8) 



where 7^ takes the big (red) loop in Fig(2J Notice that all 
the Ty k with the paths nomotopic with 7^ are equivalent. 
Also it is easy to check that T 7fe commutes with both H 
and P a . So obviously T lk \g) is another ground state of 
H. Noting that T 2 fc = 1, we can construct all the ground 
states 



i^=n^f^i°>- 



(9) 



Here, the product contains all the K "basic loops" (i.e., 
loops around each single hole). 

Note that in the representation above, T~ h \G) = ±|G). 
So we can interpret TL as the operator measuring the 
"spin" direction of the fcth hole. We will utilize this phys- 
ical meaning of T lk for encoding in the next section. 

Now we find the whole ground space of the array sys- 
tem with K holes. The dimension of the ground space is 
just 2 K . It is interesting that the degeneracy correlates 
with the topology of the system. This is why it is called 
topological degeneracy. We can encode K qubit in this 
space and do computations in it. 

Based on the knowledge of the ground space, let us 
study the low-lying excitations of the system. An open 
string product of cr z , which is compatible with the con- 
straint, excites two "charges" at the two ends of the string. 
With the help of the anti-commutation relation with the 
Hamiltonian, after a simple algebra computation, we can 
deduce that the energy of each single "charge" is 2r above 
the ground. 

Unfortunately, a string product of <j x is not allowed by 
the constraint, which is different from the case in toric 
code modelUj. To get another type of excitation, we 
should break the constraint Eq.([S} and contemplate the 
problem from the original Lagrangian of the array. Ap- 
proximately, we only keep the degree of freedom of the su- 
perconductor islands at the two ends of the rhombi (i.e., 
at the vertices of the hexagons). Correspondingly, we use 
the modified potential Eq.©. Then the Lagrangian of 



the array is 
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(^-^) 2 + 3£j cos 



A^ 



Ej cos 



3A^ 



(10) 

where A^y = ipi ~ ~ V- A ' I- Taking the continuum 
limit, the Josephson array system can be mapped to a 
Type-II superconductor [28]. A vortex solution exists in 
the system, taking one flux quanta. The energy of a 
single vortex is 
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Here, A 2 = < 5>§/97r 2 /io£'j, and £ is the length between the 
two centers of the hexagons. 

Vortices stay at the hexagons where the constraint is 
broken. The dynamics there are changed correspond- 
ingly. There is one kind of dynamic that interests us 
most, the single flip of a rhombus. To get this kind of 
dynamics, we carefully change the flux of two adjacent 
hexagrams to integer multiples of $o- Two vortices ap- 
pear and stay here. Single flip happens only at the rhom- 
bus between the two hexagons. Since the other rhombi 
stay close to the hexagons that satisfy the constraint, the 
single flip is still not allowed. This is also the reason why 
we must change the flux of two adjacent hexagons, but 
not one or two separated hexagons. If we change the flux 
of only one hexagon, we still cannot get the single flip 
because it is forbidden by the adjacent hexagon. 

Now we can add a <j x term to Hq with the amplitude 
t at the rhombus between the two adjacent vortices. Or 
rather, we can interpret it in such a way that two vor- 
tices with each energy E g are excited by a a x operation. 
This term is important for operations that we will dis- 
cuss later. So we have got two types of quasi-particle 
excitations in the array. 



III. PROTECTED QUBIT AND COMPUTATION 

After all the discussions about the properties of the 
array we focus, in this section, on the computation pro- 
cess. It was mentioned in Ref. [25j that motion of quasi- 
particles along topologically non-trivial paths changes 
the state. We care about how to achieve it. In this sec- 
tion, we show how to encode logical qubits in the array, 
and propose an exact method to operate qubits in lab. 

However, first we give a short review of stabilizer code 
formalism. In fact, this system can be regarded as sta- 
bilizer code. Then we show how the topological degener- 
ated ground space is protected from local noise. Logical 
qubits are encoded in the "holes". The spin directions of 
holes can be used to denote the state of the K logical 
qubits. Finally, we show how to implement the single 
qubit and multi-qubits operations. 
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A. Stabilizer code 

Stabilizer code is developed as a general theory of 
quantum error-correction code[l|. In this formalism, 
N physical qubits are used to encode K logical qubits 
(N > K). Certain types of errors can be detected and 
corrected so that the reliability of quantum computation 
is greatly improved. 

First, we define the Pauli group Gat on TV qubits as 

g n = {* fc K) ai K) bl , • • • , {°nY n i° 7 h) hN }, (12) 

where k,a{,bi are integers. Gn contains all the iV-body 
Pauli operators with coefficients ±l,±i. Note that any 
two elements in the Pauli group either commute or anti- 
commute with each other. 

A subgroup S of Gn is called a stabilizer group if there 
exists a nontrivial space Vs stabilized by S, 

V S = M\ g\1>) = W), VgeS}. (13) 

It can be proved that Vs is 2^-dimensional if the group 
S has N — K independent generators. We denote the 
stabilizer group as S = (gi,92, ■•■ i9n-k)i where gi is 
one of the N — K generators of the group. The elements 
in S are called stabilizer operators, or just stabilizers. 

K logical qubits are encoded in the stabilized space Vs- 
To correct errors, we should first notice that the effects of 
errors can be represented by the one-body elements in the 
Pauli group Gn, because errors are usually local. When 
no error happens, after measuring all the stabilizers in 
S (just measuring all the N — K generators is ok.), all 
results are 1. Assume one error, denoted by E € Gn — S, 
happens, so that the original computation state \tp) G Vs 
is changed to E\ip). HE anti-commutes with at least 
one stabilizer in S, e.g., g, the measurement of g gives 
the result —1, then we know the error type and we can 
correct it by corresponding operation. However, there 
may exist some elements in Gn — S that also commute 
with all the elements in Gn- This type of errors can not 
be detected or corrected. 

The basis of logical qubits should be written out for 
computation. We can always find K operators Zi from 
Gn — S such that g\ , . . . , gN-K , Z\ , . . . , Zk form an inde- 
pendent operator set. Zi plays the role of a logical oper- 
ator o~ z . The logical computational basis \xi,x% . .. Xk)l 
is actually the state stabilized by the group 

{gx, . . . , Qn—K: (-l) Xl Z 1: {-lf K Z K ). (14) 

The logical Xi can be also found from Gn — S satisfying 

1. [X u Xj] = and 

2. {Xi, Zj} = 25ij. 

We emphasize here that the logical operator set {Zi} is 
not unique. Also, we did not limit Zi to be of some 
certain form such as a product of several a z . It can be 



a mixed product of a z and o~ x , or even a product that 
contains only o~ x . The only condition required is that 
gx, - - ■ , gN-K, Z\, . . . , Zk form an independent operator 
set. For example, suppose we find a set {Zi} that sat- 
isfies the condition and {Xi} is the set of correspond- 
ing logical flipping operators, it is easy to check that 
{gx, . . . , gN-K, Xi, . . . , Xk} is also an independent set. 
So the status of logical sets {Zi} and {Xi} are the same, 
and their roles can be exchanged. 



B. Logical qubit encoding 

In this section, we show how to encode logical qubit in 
the topological protected array. As an example, the toric 
code is a kind of stabilizer code. The computation is done 
in the ground space of Htc — ~9 Yl II A s — h \~\ Bp, 
where A s and B p are the star and plaquette operators. 
The ground space is stabilized by the operators {A s ,B p }. 
The loop products of a x and a z along two directions give 
rise to the logical gate Z\, Z<z, X\ an d AV 

But toric code model only encodes two logical qubits. 
If we want to encode more qubits, we have to implement 
a much too complex boundary condition for more genus. 
In our system with puncture topology, we can see that 
encoding is much easier. 

We can write down the stabilizers in our considering 
model, {Pa} and {Q a bc}- The physical ground space sat- 
isfying P a \ijj) = Qabc\ii} — constitutes the protected 
code space. But this is not enough. Now we have to find 
the logical operators {Zi} and {-Vi}. In an array with K 
holes, we can see that the product of a x around the fcth 
hole (i.e., T lk mentioned earlier), actually plays the role 
of logical operator Zj. . We have mentioned that eigenval- 
ues ±1 of Zk can be regarded as the spin values of the 
fcth hole. We can also find the logical operator Xk as the 
product of a z along the path from the outer boundary 
to the inner of the /cth hole (The yellow lines connecting 
the right hole to the outer boundary in Fig(5]). Xk flips 
the spin of the fcth hole to the opposite direction. 

It may feel strange that we take the loop product of 
o~ x as logical Z while taking the string product of a z as 
logical X. We should notice that the string product of 
a x is not allowed by constraints at the two ends of the 
string. On the other hand, although we can exchange the 
roles of X and Z as mentioned earlier, the representation 
we chose here has a better physical meaning as "spin". 

Therefore, we denote the logical computation basis by 
the spin directions of the holes, like | 1HJJ-2 • • • 1H • ■ •)• 
The corresponding single qubit operations are Zk and 
Xk, formally. We emphasize again that all the nomotopic 
paths are equivalent. 

Now we show why the code space is protected against 
local noise. We assume that the noise can be treated 
as local perturbations on the rhombus, a 2 .. To consider 
the perturbation effect of the first order, we compute the 
matrix element (G a \o-^.\G p) . {\G a ), \Gp), . . .} is a set of 
orthogonal basis of the ground space Vg, as discussed 
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of perturbation appears in the (AT — l)th order. N is 
the number of rhombi along the path. One can get the 
energy splitting [21] 



5e = J2(1\k \H'( 



1 



Eq — Hq 



H'y\u) 



(=^I> (16) 



Figure 3: A loop product of a x around the fcth hole measures 
the direction of its flux. The state shown is | iYxlT2-U-34J'4 ) - 



in the last section. There must exist at least one Q a bc 
that anti-commutes with cr|, so aZ\Ga) is driven outside 



Vg and (G a \af \Gp) is zero. Only perturbations of an 
order higher than N take effect, such that N Pauli oper- 
ators form a topologically nontrivial path across different 
boundaries, or form a loop embedding at least one hole. 
When N is large enough, the system can be well pro- 
tected and this is the reason why it is called topologically 
protected. 

When it is cooled to a low enough temperature, the 
system arrives at the protected ground space automat- 
ically. The computation is well protected from local 
noises. However, proper perturbations can drive quasi- 
particles tunnel through nontrivial paths, and this can 
be used to construct quantum gates. We will do compu- 
tations by this method in the following part. 



C. Single qubit operation 



where Eq is the ground energy of Hq. Each H' con- 
tributes one single a z . Only the term that contains a 
complete string product of a z connecting both the inner 
and outer boundaries takes effect as X/. (remember that 
in the representation we use here, products of a z behave 
as X not Z\), while others vanish. So we get the effective 
Hamiltonian of the array as H^ ff = SeXk ■ 

In this process, a charge tunnels across the array from 
the fcth hole to the outer boundary. Degenerated states 
split just as in the case in the double- well of a single 
rhombus. 

Similarly, we want to get an effective Hamiltonian like 
^eff ~ by adding perturbations of the form such as 

& x ■ However, as mentioned before, this is not allowed 
by the constraint under ideal condition. Something more 
should be done. We have talked about the vortex exci- 
tation in the last section. To achieve this purpose, we 
carefully change the flux of the hexagons along the loop 
around the fcth hole (as shown in Fig|2]) to be integer mul- 
tiples of $o- The constraint is broken, and the dynamics 
of the corresponding hexagon changes. The single flip of 
a rhombus is allowed. So we get the desired perturbation 



We showed how the computation basis for logical 
qubits are represented. Now let us see how to imple- 
ment quantum gate operations in this system. From the 
discussion of the system previously, we know that quasi- 
particles are excited by Pauli operators. The product 
operators of a x and a z along topologically non-trivial 
paths change the state of the system. So we get a more 
physical picture that the transmission of both charge and 
vortex changes the state. This can be used for operation 
in this section. 

To complete the qubi t g ates operation, we use a 
method similar to Kou's[2f(. We change the external 
magnetic field. Effectively, the ground space behaves as 
a pseudo-spin system. We can see that it just drives the 
quasi-particles to tunnel through different sectors. 

As an example, we slightly tune the flux of the rhombi 
along a path from the fcth hole to the outer boundary. Ac- 
cording to the previous discussion, perturbations of the 
<j z form appear in the single-rhombus Hamiltonian, and 
this is allowed by the constraint. So we get the Hamilto- 
nian 



path 



H' 



(15) 



Confined in the degenerate ground space, the main 
part gives identity with an energy coefficient. The effect 



loop 



(17) 



where i takes the value in Eq.Q. Two vortices emerge. 
One is transmitted around the fcth hole, and then an- 
nihilates the other one. This process gives the Z term. 
Similar to Eq. (|Trj]) , we get the energy splitting 



6E = J2(it k \H'(-^-H') 
&o — Ho 



t 



N 



(18) 

So we have got Ti^ = SEZk- 

As we known, any single-qubit rotation can be written 

as 



U(9,(f>,j) = e -iiZ e -i<t>X e -i6z_ 



(19) 



We have got both = SeZk and 7-^ ff = SEXk, then 
we can get any single qubit gate, at least in principle [2^]. 

We can tune the corresponding flux of rhombi or 
hexagons and get the effective Hamiltonian we need, care- 
fully control the evolution time of the system (that is to 
tune the parameters 7, 0, 9 above), then we can construct 
any single qubit gate we want. 
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Figure 4: (Color online) Multi-qubit operation is sketched 
here. The (yellow) solid line represents the rhombi with flux 
tuned by 5$. This gives the effective Hamiltonian T-L%ff = 
aXm®X n . The red dashed loop represents the hexagons with 
flux changed to integer multiples of $o- That gives Heft = 

/3Z m ® Z n . 



D. Multi-qubit operation and CNOT gate 

In the last part, we saw that any single qubit gate can 
be constructed. In the this section, we show that in this 
array multi-qubit operations are also easy to realize on 
any two qubits. 

When we change the flux of the rhombi on the line con- 
necting two holes (see the yellow line in FigH]), similarly, 
the states of the corresponding two qubits are flipped 
at the same time, and this gives X m ® X n term in the 
Hamiltonian. By using the perturbation approach, the 
amplitude is 



" = £<1UU \H'( 

3 



1 



Eq — Ho 



(AE/2) N 

(_ 4r )AT-l 



Here, H' 



AE 



. (20) 

^2 o~ z , and the summation contains all 
the rhombi on yellow line. We get a effective Hamiltonian 

Analogously, we can get another type of operation. We 
can tune the flux of the hexagons along the loop around 
two holes to integer multiples of $q. This gives a pertur- 
bation H' = tJ2< jX ■ The amplitude is 



= I>m1N, \H'( 



1 



Eo-% 



t 



(21) 

Then we get the effective Hamiltonian H^ff — f3Z m (g) Z n . 

The most important thing we are concerned with is to 
construct CNOT operation, because CNOT gate together 
with arbitrary single qubit rotation gate gives rise to a 
universal quantum computation. This is also realizable 
in this array. 



First, we should construct another effective Hamilto- 
nian. We tune the flux in such a way that a vortex is 
transmitted around the mth hole, and a charge is trans- 
mitted from the nth hole to the outer boundary at the 
same time. The effective Hamiltonian contains contribu- 
tions of three parts, H^g = 71 Z m + 72 X 

n ~t~ 73 Z m <Sl X n . 

71, 72 and 73 = 7172 are evaluated in the same way with 
5E and <5e in the last section. 

A CNOT gate can be constructed with the help of 

^eff> ^eff and ^eff We Can let ^eff eV ° lve for *1 = 4^ ' 

Then we should change the external magnetic field to 



3tt 
471 



■tl 



construct ~H. z ef( and let the system evolve for t% 

Finally, we construct ji^ff and ^ the system evolve for 
<3 = — ti . Now we get the CNOT operation with an 

external global phase e~'" r / 4 , which is negligible. 



cnot — e 



'ti e -iU*t2 e -iH x t a 



(22) 



IV. SUMMARY 

In summary, we propose a TQC scheme with protected 
qubits in a JJA system, and we indicate how to perform 
qubit encoding, single qubit operation, and multi-qubit 
operation. Especially, we show the way to implement the 
CNOT quantum gate operation which is a key element 
in quantum computation. 

The scheme we propose here can be realized in real JJA 
system since the technology of the Josephson junction is 
developing rapidly. Considering the convenience of re- 
alizability and scalability of our system, the dimension 
of ground space is related to the number of punctures 
on the array. The logic qubit encoding here is in a pro- 
tected space so that our system is topologically protected 
from noise. Besides, specific quasiparticles are driven by 
proper perturbation transmitted along topologically non- 
trivial path. This makes small splittings happen. We give 
the effective Hamiltonian under different condition cor- 
responding to various splitting situations and this can be 
used to do quantum computation. 

The work is supported in part by the NSF of China 
Grant Nos. 90503009 and 10775116 and 973-Program 
Grant No. 2005CB724508. 
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